home *** CD-ROM | disk | FTP | other *** search
/ Language/OS - Multiplatform Resource Library / LANGUAGE OS.iso / cpp_libs / rwvector.lha / RWVector2.1 / src / mathpack / dcffti.f < prev    next >
Text File  |  1989-08-14  |  2KB  |  49 lines
  1.       subroutine dcffti (n,wsave)
  2. c
  3. c     Double precision version.  -tk
  4. c
  5. C***BEGIN PROLOGUE  CFFTI
  6. C***REVISION DATE  811015   (YYMMDD)
  7. C***CATEGORY NO.  D6
  8. C***KEYWORDS FFT,FAST FOURIER TRANSFORM,FOURIER TRANSFORM,COMPLEX
  9. C***DATE WRITTEN  FEBRUARY 1978
  10. C***AUTHOR  SWARZTRAUBER P.N. (NCAR)
  11. C***PURPOSE
  12. C   INITIALIZE FOR CFFTF AND CFFTB
  13. C***DESCRIPTION
  14. C     *****************************************************************
  15. C
  16. C     SUBROUTINE CFFTI(N,WSAVE)
  17. C
  18. C     *****************************************************************
  19. C
  20. C     SUBROUTINE CFFTI INITIALIZES THE ARRAY WSAVE WHICH IS USED IN
  21. C     BOTH CFFTF AND CFFTB. THE PRIME FACTORIZATION OF N TOGETHER WITH
  22. C     A TABULATION OF THE TRIGONOMETRIC FUNCTIONS ARE COMPUTED AND
  23. C     STORED IN WSAVE.
  24. C
  25. C     INPUT PARAMETER
  26. C
  27. C     N       THE LENGTH OF THE SEQUENCE TO BE TRANSFORMED
  28. C
  29. C     OUTPUT PARAMETER
  30. C
  31. C     WSAVE   A WORK ARRAY WHICH MUST BE DIMENSIONED AT LEAST 4*N+15
  32. C             THE SAME WORK ARRAY CAN BE USED FOR BOTH CFFTF AND CFFTB
  33. C             AS LONG AS N REMAINS UNCHANGED. DIFFERENT WSAVE ARRAYS
  34. C             ARE REQUIRED FOR DIFFERENT VALUES OF N. THE CONTENTS OF
  35. C             WSAVE MUST NOT BE CHANGED BETWEEN CALLS OF CFFTF OR CFFTB
  36. C
  37. C***REFERENCES
  38. C***ROUTINES CALLED  CFFTI1
  39. C***END PROLOGUE  CFFTI
  40.       implicit double precision (a-h,o-z)
  41.       dimension       wsave(1)
  42. C***FIRST EXECUTABLE STATEMENT  CFFTI
  43.       if (n .eq. 1) return
  44.       iw1 = n+n+1
  45.       iw2 = iw1+n+n
  46.       call dcffti1 (n,wsave(iw1),wsave(iw2))
  47.       return
  48.       end
  49.